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ABSTRACT 



Real time simulation of the time-dependent behavior of low power 
thermal reactors is amenable to analog computer techniques. This paper 
presents the basic physical concepts relative to thermal reactor operation, 
a derivation of the governing transient equations noting simplifying 
assumptions, and the development of the analogous electrical circuit to 

t 

ultimately provide solutions to these equations. In addition, many special 
problems involved in the simulation of a specific reactor are discussed 
and analyzed in detail. As a result, the successful design and con- 
struction of a simulator whose time-dependent behavior closely approxi- 
mates that of an actual reactor is demonstrated. The reactor simulated 
is the Aero jet -General Nucleonics low power thermal reactor (AGN-201) 
currently installed at the United States Naval Postgraduate School, 
Monterey, California. Flux-time recorder plots for various sequences of 
operation of reactor and simulator are included to indicate the degree of 
simulation attained. 
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A LOW POWER REACTOR SIMULATOR 



1 . Introduction . 

Detailed analysis of the transient behavior of nuclear reactors 
has been greatly facilitated by the use of electrical or electronic 
simulators. A reactor simulator is essentially an analog computer with 
circuitry specifically designed to solve, on a true time base, the 
rather complex differential equations governing time-dependent reactor 
operation. 

As will be evident from the development of the governing equations, 
complete solution by direct analysis would be extremely laborious, if 
not impossible, for other than a few limited cases where modifying 
assumptions must be made before a useful solution can be obtained. 
Moreover, direct mathematical analysis could not be expected to instill 
in one any real feeling for reactor kinetics under typical operating 
conditions. Because circuit values can be measured quite precisely, 
accurate and continuous solutions to the transient equations for almost 
any imposed operating conditions may be readily obtained from the sim- 
ulator. In addition, by employing a simple simulator control system, 
a strong feeling for reactor responses to any of a multitude of typical 
operating situations may be achieved. This particular feature of the 
simulator might permit reactor operating personnel or students to receive 
the equivalent of a great deal of reactor control experience and 
operational training when access to an actual reactor is not feasible . 

Because analog computers are well-tried and accepted engineering 
tools for use by those involved in the study and design of dynamic systems 
of all types , reactor simulators find application in preliminary reactor 



1 



design and especially in the evaluation and perfection of associated 
systems which might be ultimately employed to control an actual reactor. 

In other words , completed control components and systems may be connected 
to the simulator by means of electronic transducers thereby permitting 
operational checks to be made on their specific design functions prior 
to actual installation. 

Depending on the complexity of the system to be studied and the 

degree of simulation desired , simulators may range in size from small 

units such as the one discussed herein to extremely large, complex, and 

expensive systems involving the use of great numbers of amplifiers, function 

generators, servo and electronic multipliers, electro-mechanical devices, 

and countless precision components. In the case at hand of a particular 

low power thermal reactor, utilization of the basic simulator circuit 
% 

together with a relatively simple control network simulates actual reactor 
operation to a surprisingly good degree of accuracy. 

The discussion of basic concepts , the derivation of governing 
equations and the build-up of the electrical analog represent a compilation 
of treatments given on this subject by the first seven authors listed 
in the bibliography. The simulator itself is essentially a Pagels £5^ 
simulator, modified to duplicate more closely the responses of the 
individual reactor under study. 

It must be remembered that the degree of simulation attained can be 
no better than the degree to which the differential equations utilized 
define the system under study. The transient equations employed in 
thermal reactor simulation contain many simplifying assumptions and 

^Numbers in brackets refer to Bibliography in Section 8 . 
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approximations} and as a consequence, precise duplication of actual 
reactor behavior nay not be expected. Despite these assumptions, however, 
it will be seen that simulator solutions to the equations employed do 
closely approximate observed reactor behavior. 

It is the purpose of this thesis, therefore, to present a rather 
complete study in the design of a reactor simulator with the ultimate 
goal being to construct such a device and then compare time dependent 
behavior of the actual with that of the analogous system. 

2, Basic Thermal Reactor Concepts. 

In order that the discussion of the simulation of a low power 
thermal reactor be complete, it is important that basic concepts be 
briefly reviewed prior to the derivation of the governing differential 
equations . 

When the nucleus of a fuel atom fissions, it splits into two primary 
fission fragments of more or less equal mass accompanied by the prompt 
emission of from two to three neutrons of high energy. On the average, 
about 2.5 neutrons are liberated per thermal neutron absorbed in fission 

C 

of a U 23 nucleus. The majority of these have energies of about one or 

two Msv. Most of the fission fragments thus formed are radioactive and 

some are themselves high energy neutron emitters. The decay constants 

for six delayed neutron precursors have been measured although the precise 

87 89 

origin of only three of the six is currently known, namely, BR , BR , and 
I 137 . The average energy of the delayed neutrons is about .5 Msv, The 
properties of the delayed neutrons play a vital role in time-dependent 
reactor behavior and permit safe control of nuclear reactor fission 
processes which might otherwise be extremely difficult. 

In addition to having small thermal neutron absorption cross-sections, 
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good moderator materials are those made up of elements of low mass number, 
where the energy transfer from neutron to moderator nucleus on collision 
may approach 100^ as in the case of a hydrogenous moderating medium. 

This results in extremely rapid slowing down of the high energy neutrons 
to so-called thermal energies (about ,025 ev) , At this point, the neutrons 
are essentially in thermal equilibrium with the moderator atoms (or mole- 
cules) and hence are referred to as thermal neutrons. Despite the 
extremely short slowing down time, some of the neutrons will either 
escape from the reactor or be absorbed by the core material before 
thermalization so as to be unavailable for subsequent fissioning or neutron 
regeneration , 

Once the fast neutrons have been reduced to thermal energies, they 
may leak or diffuse out of the reactor or be absorbed by moderator, fuel, 
or poison nucleii. A large fraction, about & 5 %, of the thermal neutrons 
absorbed by the U 23 fuel will cause fission; and a new neutron lifetime 
cycle begins again. The next job, therefore, will be to express in equation 
form, the lives of neutrons as described above. 

3. Derivation of the Transient Equations. 

Based on the conservation or balance of neutrons principle, applicable 
to any reactor, a general expression for the rate of change of neutron 

population density (n) with time can be written as follows; 

bn 

Production-Absorption-Leakage = — (1) 
The rate of production of neutrons may te obtained from empirical con- 
siderations . Likewise , the rate of absorption depends on various absorp- 
tion cross-sections which must be determined experimentally. For determin- 
ing leakage, resort is made to classical mechanics and the familiar laws 
of diffusion theory which, although somewhat more complicated, can be 
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reduced to simplified expressions when certain prudent assumptions are 
introduced . Let us now form the basis on which the various simplifying 
assumptions can be made. 

The thermal reactor to be considered here is of the homogeneous type 
with reflector. The fuel is a homogeneous mixture of U 23 ^ in the form of 20«g 
enriched U0 2 with polyethylene, (CI^ n , as the moderating material. A thin 
high density graphite reflector is employed. However, in the analysis 
which follows, it will be assumed that the reactor is homogeneous but un- 
reflected. The degree to which this assumption is valid can best be 
appraised by a comparison between final simulator solution and actual 
time-dependent reactor test runs. As will be seen later, the assumption 
was quite reasonable in this case. 

Before continuing with additional assumptions , it will be helpful to 
rewrite equation (1) as the familiar diffusion equation for thermal 
neutrons in non-equilibrium, basic to reactor theory, namely 

S - 2 a 0 - (-DV 2 <£) = ^ (2) 

where S is the production term in neutrons per cm 3 per second , 2 is the 
macroscopic absorption cross-section for monoenergetic (or thermal) neutrons 
in cm \ 0 is the. thermal neutron flux in neutrons per cm 2 per sec, 'and 
(— D V is the neutron leakage per cm 3 per sec. By the acceptance of 
equation ( 2 ) as the equation representing (l), it has been assumed that we are 
dealing only with monoenergetic (or thermal) neutrons and then only at dis- 
tances greater than two or three mean free paths from strong absorbers, 
sources, or boundaries. This also assumes that the moderating or diffusion 
medium is a poor absorber relative to its ability to scatter and, hence, slow 
down the high energy neutrons. The partial derivative of n is utilized because 
n is a function of both position and time. However, it will be presumed 
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that the reactor is at or near critical and the shape of the flux 
distribution remains constant even though it changes in amplitude with 
time. Under the latter provision, (-DV 2 #) could be represented by some 
constant multiplied by the flux <t>. feking use of this, therefore, let 

(— DV a 0) = 2 where the symbol 2 will be used throughout this 

1 

development to indicate the macroscopic cross-section for or the probabil- 
ity of occurrence of the particular process indicated by the subscript, 
i.e., in this case 2 ^ to mean probability of loss by leakage or leakage 

coefficient. This essentially eliminates the space variable and permits 

cln £)n 

the use of the total derivative of n, or in lieu of ^ The flux 4> 
under the above considerations is merely the product of n, the density 
of thermal neutrons in neutrons per cubic centimeter, and v, the mean 
thermal neutron "speed." Or, expressed in equation form, 

<P = nv neutrons /cm 2 /sec 



Thermal neutron flux therefore, is a scalar quantity. We are now in 
a position to examine each of the terms of equation (2) in an attempt to 
develop a more workable and meaningful equation. 

The source term (S) must be broken into several parts, i.e., S p , 
and Sq , referring respectively to prompt , delayed , and those uniformly 
distributed as a consequence of an extraneous outside source (e.g., a 
Ra-Be source). Expressed in equation form; 
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S = s p + S d + So (3) 

The rate of fission can be expressed as 2^0. If ~V ^ represents the total 
fission yield of fast neutrons from both delayed and prompt sources and 
P represents the total fraction of these neutrons which are delayed, then 
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the prompt neutron generation rate is given by 



2^/ (l-p) 

However, not all of these neutrons reach thermal energies. As previously 
mentioned some leak out or escape and others may be captured by fuel or 
other nuclei! during the moderating or slowing down process . Probabilities 
have been evaluated for each of these processes . The former is given by 
the fast non-leakage probability, e (where B 2 is the so-called 

buckling factor related to size and T is the Fermi Age, a measure of the 
distance neutrons travel in slowing from fission to thermal energies), and 
the latter by p, the resonance escape probability. For simplicity, these 
factors will be combined into a single thermalization probability symbolized 
by p.^ • Therefore, the prompt thermal source becomes 

S P = » 2 t tt-e) 'V' p th (4) 

-B 2 T? 

where p^ = pe . Now the remainder of the neutrons are the so- 
called delayed neutrons associated with the particular fission fragments 
mentioned in section 2 . The delayed neutron generation rate, therefore, 
is given by 

z f 0 p V 

Now, as previously stated, six delayed neutron precursors have been 
determined so that the delayed neutron generation rate of each of the six 
may be expressed by 

g A = Pi 

where p^ = fraction of total number of neutrons contained in the i 
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delayed group. It should be noted here that fs 




Pi- 



At any time , 



there will be a concentration of the i^ group of delayed neutron 

emitters present in a unit volume of the core . Because these delayed 
neutrons are emitted as a result of the normal radioactive decay of the 
precursors, each emitter obeys the familiar radioactive decay law. In 
other words , the rate of change of the concentration of the i^ group will 
be equal to the production rate minus the decay rate , or 



dt g i " Vi 

where is the decay constant for the V 1 emitter. 
Therefore , 




2^- x i c i 



( 5 ) 



Delayed neutrons are therefore supplied to the reactor at the same rate at 
which decay of the precursors takes place. This can be expressed as a 

4 

summation \^i* Because the average energies of the delayed 

neutrons are somewhat less than those of the prompt, there will be less 
tendency, for the delayed neutrons to leak out during the slowing down 
process | and hence they will have larger values of non-leakage probabilities. 
The difference in the' values of non-leakage probability for the two types 
may be taken into account by utilizing effective values of in the 

final equations which are slightly larger than the actual values, This, 
then, permits multiplication of the delayed neutron term by the same value 
of non-leakage probability as that used for the prompt neutron term. 

Combining this factor again with the resonance escape probability p, 
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the rate at which delayed neutrons are supplied can now be expressed as 
a summation of terms , namely 

6 

S d = p th \ C i ^ 

i=l 

where XjC^ is defined by equation (5). Finally, Sq comes from an 
essentially constant extraneous source of neutrons which makes a contribution 
to production rate in the form of a uniformly distributed source . This 
term will be discussed in more detail later. 

Ihe absorption term in equation (2) needs no additional explanation 
except to say that Z & represents the macroscopic absorption cross- 
section for all processes , that is absorption of thermal neutrons by 
structure, moderator, fuel, or poison nuclei!. 

The leakage term was discussed earlier and can be taken to equal 
where is dependent on the moderator characteristic and 

reactor size . 

The component terms of equation (2) can now be assembled as follows; 



dn 

dt 



(s ♦ S d ♦ s 0 ) - z a * - Z 1 * 



dn _ 
dt 



* 2 f(l - P )V Pth * *th * s 0 

i=l 







or 



ft * i 1 -¥ * Jj V^] - *( 2 . * 2 ^ * So (7) 



In order to reduce equations (5) and (7) to more usable forms, several 
terms must be defined. First, the term reactivity jO may be defined as 




Jk£± 



K 



eff 



( 8 ) 
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where is the effective multiplication factor for a given finite sized 

reactor. By effective multiplication factor is meant the population ratio 

£ 

of successive neutron generations. If 1 is next defined as the mean lifetime 

* 

of thermal neutrons in the reactor and lA-o as the rate at which neutrons, 
prompt and delayed, are thermalized in the reactor, then could also 

be expressed as the ratio of 1* /l^ 56 ( , or 



K eff 



i A 0 



(9) 



This may be interpreted as follows: is the ratio of the time spent 

by the thermal neutrons diffusing in the reactor to the time spent in being 
thermalized. In other words, if it takes longer to thermalize the second 
generation of thermal neutrons than the life expectancy of the first 
generation, the chain can be expected to die out. In this case, would 

be less than unity. The reverse would be true if were greater than 

unity . When K eff is precisely unity (or in the case of the presence of 
an extraneous source, slightly less them unity), just as many thermal 
neutrons are born as die and the chain is sustained. These three situations 
represent subcritical, supercritical, and critical behavior of the reactor. 
1^ cem be expressed as the reciprocal of the thermalization rate f>3 , or 





This term is also given the name of modified lifetime by Murray ^3 j . 
Now thermal neutron lifetime in an infinite medium would be given by the 
average distance traveled by the neutron before absorption, called the 

absorption mean free path, X ( = — ), divided by the mean neutron 

« a 

velocity v. However, in a finite reactor, this must be reduced by the 



( 10 ) 
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thermal leakage factor 



Thermal Absorption 



Thermal Absorption + Thermal Leakage 



or 



a 



vZ 2. + Si " ,( 2. *Zi> 



(n) 



Substituting (9)* (10), and (11) into (P), therefore, 



/° 



= 1 - 



l*Al 



= 1 - 



v( 2 a + JZj) 

p th 



( 12 ) 



Now, utilizing 1*, l£ and as defined above, equations (5) and (7) may- 
be simplified as follows; 



& - U - f>K * Z Vi>pJ| - * < 2 a ^ 



dn 

dt 



) + s c 






(1 - p) + P 



th 






♦ S r 



■f” ( 1 - P - ^ + P+h \C. + s c 

it eff th i=i i i 



0r dt P n + p th X i°i + S ° 

J.Q I-± 



(13) 



^ dC. 

Likewise, from (5), ^ = ( Z f ^ ) 



or 



Pi 



dC. 



lo p th 



n - 



dt 



( 14 ) 



Equations (13) and (14) are the desired transient equations. 
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Finally, a single governing equation can be written by a simple 
substitution of (14) into (13) as follows j 



It will be further assumed that the so-called one-group theory is 
sufficiently accurate for accomplishing the goal at hand. In one group 
theory it is supposed that all production, diffusion, and absorption of 
neutrons occur at a single energy, namely, the thermal energy. Although 
this is a considerable approximation, the final simulator circuit is 
thereby simplified and results are still quite satisfactory. Since this 
theory implies that all fission neutrons have thermal energies, the 

resonance escape probability, p, and the non-leakage probability during 

-B 2 T 

the slowing down process, e , do not affect the multiplication factor 

and, hence, become unity in equations (13) thru (15). The two' terms are 
carried through in the development of the electrical analog # but no 
attempt is made to employ them as variables in the final circuit. 

4. The Electrical Analog. 

Equation (15) is now in a convenient form for a computer solution. 

The network given in Figure 1 will be shown to describe this differential 







1=1 



(15) 
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FIG. 1. BASIC REACTOR SIMULATOR CKT. 



1 *3 



equation. Referring to Figure 1, and using Kirchoff's second law, it 
is evident that 



h = 




+ 



i 
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Note the similarity in form to (15). 



The Rate Term 



( 1 * — ) 
(1 ° dt > 



= c 



dE 
7 dt 



l 



and letting E in volts be proportional to n in neutrons per cm 3 per 
sec by the constant 3 , i.e., E = 3 n, then 



I 
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dn 

dt 



(16) 



The Reactivity Term ^71) 

Let A be the proportional distance of the potentiometer tap (at P) 
from the center of the potentiometer. Referring to the more detailed 

t 

development given in Appendix I, 



R |I + R/P + P/4R (1 - A 2 Fj 

The range of values of A are zero at the center, + 1 at the + E end, 
and -1 at the -E end. Now, it is apparent that Ig will have its 
maximum positive value at 4 : + 1, Therefore, the range of reactivity 
variation of the simulator may be set and the proper scaling factor 
determined by setting 



14 



s- 



f' 



max 



8 

8 max 



T 



or 



Z 5 



1 * 4R(R ♦ p) ^ 






max 



(17) 



Since Ig occurs at A = 1 » then 



max 



E 



8 max 



R(l + R/P) 



and therefore, 



R(1 + R/P) 



t4 



r 



max 



igain representing E by vf n, this reduces to 



8 



J 

r(i + r/p) 



—}/° n 

max*-'/ 



(18) 



Nfore careful consideration of (17) above reveals that jO is a non-linear 
function of A which is undesirable for our purposes . This is due to 
the potentiometer loading having its greatest effect around the mid-position 
of the potentiometer. Linearity can be very closely approximated, how- 
ever, by making R sufficiently large in relation to P that the 
denominator of (17) becomes essentially unity, leaving 



MX 
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If, for example, R=5P» the linearity remains below one percent for any 
value of A from (-)l to (+)1. 




The Delayed Precursor Term i=l 



. v 

The term of (15) is simulated by the six RjC^ delay networks. 

The basic voltage equation describing the voltage drop across any of the 
i^* 1 type networks is 



E = 





where is the capacitor charge in coulombs and c^ is the capacitance 

of the i^ 1 group. Ihe electrical charge is analogous to the "storing" 
of delayed neutrons in the fission fragments before they are emitted. The 
corresponding delayed neutron equation represented by these networks is 
(14) , which upon rearranging becomes 



n = 



p th X 1 C 1 

Pi 



^ Pth 

dt 



Pi 



Again letting E = n and matching terms with the basic equation above , 



^5 1 0 P+w i 

Qi = ~h Cl and R i c i = ^ 



Now, I. = 



dQ a 

IT 



( A \ i* dc. 

( n D ) l0 P-th * 



M 



i“i 



dt 



dC. 



and by substitution for from (14) , 
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^ - ( ' ?1 " " ( ^ ^ Pth Vi 



Pl R l 



But, I, = If + I 2 + I 3 . . 

1 



and 



l ?a] p th **- = Pi “ " [lA - ] ^ f5 



x i c i (19) 



Therefore 



» “ P» ~ P2 R 2 = ?3 R 3 



The Source Term (S.) 



I f > > P 9 , then 



6 9 V 



"9 



where A^ is the fraction off ground potential and V is a constant 
voltage source. Now, letting A^ = 



-So- 



So 



(max) 



then 



*9 = 



*9 S ° (max) 1 * 



i* S 
x o ^0 



( 20 ) 



Finally, summing the currents as given by equations (16), (18), (19) 
and (20), 



X 7 = *8 



O 

-2 It 



+ Ir 



i=l 
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[J t 





In order for the current output of each term to be of the correct 
magnitude, the scaling factors which are enclosed by brackets in 
equation (21) must be identically equal to each other; that is, 





J _ 




( 22 ) 



If this can be done, then these coefficients drop out of equation (21) 
leaving equation (15) which is to be duplicated electrically. Values to 
insert into (22) will be determined by constants in the actual system to 
be simulated as well as arbitrary values selected in accordance with standard 
computer operating techniques. 

5, Simulation of a Specific Low Power Thermal Reactor — AGN-201 

General thermal reactor simulation, therefore, can be readily duplicated 
by the circuit given in Figure 1. This circuit when set up on a conven- 
tional analog computer, using completely arbitrary circuit parameters, will 
yield typical time-flux curves. The problem at hand, however, is to 
simulate a spec ificre ac tor wherein lie many additional considerations which 
first must be resolved before a final circuit can be drawn up. These 
include detailed investigations to. 

(1) determine actual range of reactivity variation, compute rates 
of reactivity change for control rod movement, and provide suitable means 
for duplicating control rod action; 
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(2) determine a means of introducing a maximum reactivity control 



and temperature coefficient of reactivity; 

(3) solve electrical problems resulting from real time simulation over 
wide operating ranges of neutron flux from source level to full power; 

(4) and to determine an appropriate value of neutron source voltage 
corresponding to the extraneous source production rate. 

In the analysis of these problems, attention shall be restricted to 
the simulation of the low power thermal reactor currently installed at the 
United States Naval Postgraduate School. This reactor will hereinafter be 
referred to as the AGN-201. 

Before solving the problems, it is necessary to mention a few essential 
features of the AGN-201 which will form the basis on which ultimate 
simulator construction may proceed. The general characteristics of the 201 
reactor are listed in Table 1 below. 



TABLE 1. GENERAL CHARACTERISTICS AGN-201 



Power (licensed ) 



100 miliwatts 
4.5 x lO^cra*" 2 sec”^ 

665 gm U 23 '* 

10" dia. x 9^" high 
11 kgm 

20 cm thick (1.70 gm/cm 3 ) 



Nfaximum Thermal Flux @ 100 mw 



Critical Mass 



Core Size (cylinder) 
Moderator (polyethylene ) 
Reflector (graphite) 
Shield (lead) 



10 cm thick 



(water) 



55 cm thick (1000 gallons) 
6£ ft 
9 ft 

-2.75 x 10 -4 /°C 



Over-all diameter 



Over-all height 



Temperature coefficient of reactivity 
(based on water shield temp, data) 
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Figure 2 shows the AGN-201 core tank containing the uranium-polyethylene 
core and portions of the graphite reflector. Figure 3 is a cross-sectional 
view of the entire reactor unit showing the core, reflector, lead and 
water shields, glory hole, and control rod assemblies. The critical mass 
consists of about 6f5 grams of U 23 in the form of 20$- enriched U0 2 , It 
is contained in the form of a powder dispersed in a solid polyethylene 
moderator. The core is made in nine sections of various thicknesses 
(Figure 2). Each of the four bottom discs has three 5-cm holes and one 
2.5-cm hole. The 201 has two safety rods and two control rods. Ihe 
three larger holes are for two safety rods and one coarse control rod while 
the smaller hole is for a fine control rod. The two safety rods and the 
coarse control rod are identical and each contain about 14 grams of U 23 . 
The fine rod containsabout 3.5 grams of U 23 "\ All rods operate so that 
reactivity is increased as the rod is inserted . The amount of reactivity 
each rod controls increases as the amount of active material contained 
increases. With the above amounts of U 23 "* in the rods, the safety and 
coarse control rods control about 1.1$ reactivity while the fine control 
rod controls about 0 .% reactivity. The rods are inserted by a reversible 
d-c motor drive mechanism requiring about one minute to travel the full 
25-cm stroke. Referring to Table 1, the thermal flux at licensed full 
power is given as 4.5 x 10^ neutrons /cm 2 /sec. This is well below that 
at which Xenon poisoning becomes a consideration, namely, at a flux in 
excess of about 10 neutrons /cm 2 /sec. Finally, because of the low 
power developed, heating of the core is so slight that it is essentially 
ambient temperature which determines core temperature and, hence, the 
maximum reactivity available . The negative temperature coefficient listed 
in Table 1, therefore, is based on the change in temperature of the 1000 
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gallon water shield surrounding the core tank £loJ • With these features 
in mind, let us proceed with the investigations into the special problems 
mentioned . 

First, a control rod calibration must be made , the result of which is ^ 
shown graphically in Figure 4. The term relative reactivity as used here 
implies an arbitrary zero reactivity reference when the rods are fully 
retracted. As will be seen later, the major portion of the reactivity 
contribution or "worth" of the rods is employed in reducing the amount 
by which the reactor is subcritical. Only a small fraction of this 
relative reactivity is actually positive reactivity as defined by (8), 

This curve represents a superposition of data taken on both fine and 
coarse control rods which actually exist on the reactor, thus providing 
a single curve to represent the movement of both rods simultaneously. 

In view of the fact that the action of the two control rods on the reactor 
can be duplicated more simply and more economically on the simulator by 
a single controlling device , simulation of operation of each rod individ- 
ually does not seem warranted. Now, under normal conditions, the reactor 
will be operating in the straight line portion of the curve when near- 
critical conditions exist. As indicated in Figure 4, the slope of this 
straight line is 837 jy© (yOx 10~^) per cm of rod travel. From weight 
data on the fuel discs which make up the two safety rods and the fine 
and coarse control rods, it is found that the relative reactivity con- 
tribution of each safety rod is essentially equal to that of the coarse 
control rod or about 11,275 p/o each. Therefore, the total contribution 
in terms of relative reactivity of all rods is found to be 36,725 . 

The reactivity potentiometer can now be designed to match the control and 



safety rod curves. Figure 5 shows a plot of percent of total resistance 



of the reactivity potentiometer versus percent of total rotation (about 
350°) where all rods are being considered. Notice that the curve is 
shown up to the point at which the safety rods are completely inserted 
as a dashed line. Because the safety rods are not utilized for control 
purposes, their individual calibration curves are of no concern as long 
as their ultimate contribution to total relative reactivity is known, 
namely, 22,550 or about 61.4$ of the total available for addition 

by rod insertion. From this point on, however, it would be desirable to 
follow the remainder of the curve which is essentially Figure 4 super- 
imposed but plotted using non-dimensional relationships. To duplicate the 
entire curve precisely would require a non-linear potentiometer which could 
be designed but at considerable expense. Instead, in view of the fact that 
the straight line portion of the combined control and calibration curve 
is of greatest interest, a center-tap potentiometer permits precise duplica- 
tion in this region. At 50$ rotation, extrapolation of the straight line 
portion back to this point reveals that 25$ total resistance should result 
at this point. For the 20,000 ohm potentiometer selected, this requires a 
5000 ohm resistor in parallel with the first half of the potentiometer. The 
arrangement is shown schematically in Figure 5; and the actual plot of 
percent resistance versus percent rotation resulting is indicated by 
circles in the figure. It is noted that 6l.4$ of the total resistance 
will now occur at about 75$ of total potentiometer rotation. If this point 
is denoted for the simulator as having a control rod position reading of 
zero, then a discrepancy will exist between the control rod- re activity 
relationships of the actual reactor and that of the simulator. As seen 
from Figure 5, the simulator rods increase reactivity linearly over their 
entire travel while actual rod insertion increases reactivity linearly 
only over the last 10 cm of travel. This situation must be accepted when 
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the actual reactivity curve cannot be duplicated over its entire range. 
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The rate of potentiometer rotation is easily obtainable by timing the 
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actual control rod insertion rate. Both fine and coarse control rods have 
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been set to move from the fully retracted to the fully inserted positions o~ 
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(0 to 25 cm.) in about one minute at the fast speed setting. Therefore, 



the last 10 cm of actual rod travel will be completed in 0.40 min. From 
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Figure 5 the last 10. -cm of actual rod travel commences at 85$ of total 
potentiometer rotation or, alternately, over the last 0.15 revolution. 
Hence, the reactivity potentiometer drive rate should be 0.15 revs/0,40 
min = .37 rpm, at the proper rate. Over this straight line portion, 
reactivity change at fast speed is found to be 



83 7 jyo /cm x cn/sec = 349 jyo /sec 

To solve the second problem, a means of establishing a variable zero 
or critical reference level must be provided for the reactivity potentio- 
meter. 'Hie reactor functions between certain definite reactivity limits 
for any specified core temperature . In other words , some scheme is re- 
quired to duplicate the normal operating ranges to be expected from 
maximum subcritical conditions to maximum supercritical conditions. For 
example, at 23°C, the maxi mum positive reactivity attained with all rods 
fully inserted is only about .0020. At the same time, with all rods fully 
retracted, the reactor has a negative reactivity of about .035. This 
one-sided variation is representative of the type of reactivity behavior 
encountered under typical operating conditions. It should be noted here 
that the above values were determined with the so-called "glory hole" empty. 
The glory hole (Figure 2) is a small access tube about one inch in diameter 
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which passes through the core center, and into which additional fuel, 
moderator, reflector material, samples and the like can be inserted to 
yield changes in reactivity. As an example, by insertion of the proper 
amounts of moderator, reflector, and lead shielding in the glory-hole, 
the core is made essentially homogeneous; and an additional reactivity 

* 

contribution of about 2600 can be realized. Moreover, an average 

linear temperature coefficient of reactivity, based on water shield tem- 
perature, of about -2.75 x 10“^ per °C exists for the AGN-201. This rather 
significant effect on reactivity must also be taken into consideration. 
Figure 6 represents the scheme devised to yield the proper operating range 
of reactivity from a relatively large negative value to a small positive 
value and, in addition, to permit control over maximum positive reactivity 
available. A permanently imposed offbalance in the values of R, and R 2 
will electrically offset the mechanical mid-position as indicated. The 
values given in the figure are those based on the reactivity conditions 
at 23°C and employing a homogeneous core. The critical or zero position 
is now located at point C. In addition, it is noted that C can be further 
controlled by movement of the maximum reactivity control and temperature 
coefficient (MRCTC) potentiometer. If, for example, the temperature in- 
creases (or alternately, if a sample with a large thermal neutron absorption 
cross-section is placed in the glory hole), the MRCTC potentiometer is moved 
as indicated in the figure. This increases the resistance in loop 2 and 
decreases the resistance in loop 1, causing an unbalance in the feedback 
circuits, therety producing a decrease in reactivity corresponding to the 

amount of resistance change applied. If, therefore, it is desired to hold 

the reactivity at zero, the reactivity potentiometer must be moved in such 

a direction as to balance the resistances in the two loops. In this case, 
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it would require movement of the arms as indicated. It can readily be seen 
that if the temperature coefficient is offset far enough, the maximum 
positive reactivity available can be reduced until the reactor may no longer 
be brought critical even though all rods are inserted. This corresponds 
to the point marked B in the figure. The original electrical center position 
C, therefore, effectively walks up and down the reactivity potentiometer 
depending on the setting of the MRCTC potentiometer. Of course, maximum 
positive reactivity may likewise be increased to prompt critical conditions 
and above if desired. One of the solutions of the transient equations for a 
step change in reactivity is the so-called Inhour equation. The plot of this 



equation using the value of neutron lifetime applicable to the AGN-201 is 









given in Figure 7. Calibration of the MRCTC potentiometer can be made by 
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easuring the stable period (T) for any particular potentiometer setting 
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when all rods are fully inserted. The corresponding reactivity may then be 

r » 

picked off an appropriate plot of the Inhour equation using the neutron life- 



/ 



time applicable to the reactor being simulated. Such a calibration curve is 
given in Figure 8 for the particular MRCTC potentiometer employed here. 
Temperatures shown in the figure are based on experimental data employing a 
homogeneous core and using the negative temperature coefficient of reactivity 
from Table 1. In summary, therefore, the MRCTC potentiometer functions as 
an electrical stop on the upper limit of the reactivity potentiometer which 
provides a scheme to preset any desired maximum reactivity (or alternately, 
any desired core temperature). The effect on the electrical circuit can be 
taken into account quite simply as is indicated in Appendix II. 

The third problem can be solved by employing stablized operational 
amplifiers with virtually zero drift to give reasonably accurate and repro- 
ducible solutions at the extremely low starting or source voltages encountered. 
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FIG. 7. PLOT OF THE INHOUR EQUATION 
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The latter are of the order of magnitude of a few millivolts if scaling is 
such as to have about four decades of multiplication from source level to 
full power. The maximum output of the amplifiers is 100 volts. A variable 
biasing arrangement is also desirable to ensure zero offset at zero voltage 
input to the amplifier gird. Any one of several suitable methods recom- 
mended by the manufacturer of the amplifiers will, satisfy this requirement. 
Causing additional problems at low operating voltages are the condensers in 
the circuit. It is found that an erroneous flux-time curve results in the 
low voltage region upon the decay of the neutron flux when the reactor 
simulator is scrammed. Analysis of oil, paper, mica, electrolytic, and 
other types of commonly used capacitors reveals a marked dielectric ab- 
sorption effect by which these types hold their charge for a relatively 
long period of time even after direct shorting across the terminals. This 
problem is virtually eliminated by the use of polystyrene or nylar condensers 
which exhibit characteristically low dielectric absorption properties. 

The fourth and final problem is significant only in that it is con- 
sidered desirable to simulate reactor responses over the entire operating 
range. This demands computer solutions to the controlling equations at r 
extremely low values of neutron flux which corresponds to the condition in 
the reactor when all rods are retracted, and the existing neutron density is 
dependent on the extraneous neutron source strength. This problem may be 
approached as follows: From equation (15), namely. 




dn 

it is apparent that under steady— state conditions -rr and terms 



drop out leaving 
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0 



n + lo s 0 = 

or rearranging, 

7° n 



( 23 ) 



And from equation (8) and (9), this can be written as 



So = 



n(i - K eff ) 



(24) 



From manufacturer's data or from actual experimental flux determinations, a 
value of n can be calculated at a level where reactivity is known. The 
latter is obtainable from experimentally determined control rod calibration 
data. Employing (23) > this will establish a value of S 0 1 q at one particular 
control rod position. From the known value of reactivity at this position, 
K eff b© determined from (8). Finally, from (24), SqI* is obtained. 

If the value of estimated neutron lifetime as given by the manufacturer is 
then employed, S 0 is readily calculated. It is noted from (22) that 



z * J. 

V« (max) R(1 * ^Y max 

or rearranging , 

T _ 1 *9 (S °(max) ^ > 

v * *u— 

Therefore , to obtain the neutron source voltage corresponding to a particular 
source generation rate, the product (So( max ) lo ) need only be determined 
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as accomplished above, and actually the two terms need not be separated 
to determine the voltage , The magnitude of neutron lifetime itself is 
important in the determination of a proper value of the integrating 
capacitor c^. Prom (22) 




or, essentially, c^ is proportional to neutron lifetime. However, if 
the order of magnitude of neutron lifetime is known, this is sufficiently 
accurate to yield excellent results as long as maximum reactivity values 
are limited to .005 and under. This figure of .005 is arrived at by evalua- 
ting the stable period corresponding to various values of positive re- 
activity when different neutron lifetimes are considered. The Inhour 
Equation may be employed to make this comparison analytically | and it will 
be found that, at low reactivities, period is independent of neutron life- 
time for a specified reactivity. Furthermore, in view of the fact that the 
stable period corresponding to a neutron lifetime of 10 sec and a re- 
activity of .005 is about 2.5 seconds, this condition on the actual 
reactor cannot be duplicated because a period scram is set to function at 
about 3.0 seconds. As a consequence, comparison runs at higher values of 
reactivity to better evaluate neutron lifetime are not feasible. 

Having solved the above four problems , the remaining circuit para- 
meters may be calculated as indicated in Appendix II and the final AGN-201 
simulator can be constructed. Figure 9 is a schematic of the final analog, 
and listed in Table 2 below are the various values of circuit components 
referred to in the figure. Also listed are the other factors which appear 
in the scaling factor equation (22) . 
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FIG*. 9 . AG-N “201 ELECTRICAL ANALOG- SCHEMATIC 
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TABLE 2a. CIRCUIT VALUES 



p 


— 


13.0 


U 


(actual value 


* 


= 


119 


Ka 








R* 


= 


129 










Cl 


= 


.0007 




R, 




102M 


C* 


= 


.0207 




Ra 


= 


30M 


°3 


s 


.208 


Pf 


R> 


= 


10. 6M 


°4 


S 


.544 


Pf 


R 4 


- 


12. OM 


C 5 


x 


2.06 


pf 


R 5 


3 


15. 4M 


c 6 


s 


.785 




"6 


x 


102M 


°7 


= 


.0039 










R g 


•& 


R f 


= 1M 








*9 


= 


10M 










V 




51 mv 











TABLE 2b. ADDITIONAL SCALING FACTOR VALUES 

* 23. OK XI (use for P in equation 22) 

R * 119KJ2. (use for R in equation 22) 

/^mx= .034775 ® 23°C (core homogeneous) 

p, = .00025 

Pa = .00085 

Pj * .00241 

P 4 = .00213 

p 5 = .00166 

P 6 = .00025 

= 2.22 

Sq = So = .579 neuts/cmVsec 
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1* 9 ^ 1 .0 x 10~^ sec 9 ^ 1^ near critical 



(i.e. K eff ^ 1) 



Figure 10 is a photograph of the completed simulator with linear 
recorder unit attached. Important features and actual circuitry are 
shown in Figure 11 which is a rear view of the simulator console . At 
the left center of Figure 11 is the reactivity potentiometer and drive 
motor assembly; at the upper right portion of the figure is the MRCTC 
potentiometer; and at the bottom portion of the photograph is the entire 
simulator circuit which is constructed to resemble the schematic of 
Figure 9. Along the top of the console are located the operational 
amplifiers and stabiliting amplifiers. 

6. Degree of Simulation Attained 

Figures 12 through 14 represent comparison flux-time plots taken 
from the reactor and the simulator, utilising the same linear recorder 
so as to ensure an identical time base. The sequence of events describing 
each of the comparison runs is briefly discussed in the succeeding 
paragraphs . 

The sequence of events describing Figure 12 is given in Table 3. 

TABLE 3. SEQUENCE OF EVENTS FOR FIGURE 12 

TIMS IN REACTOR SEQUENCE CORRESPONDING 

MINUTES SIMULATOR SEQUENCE 



Steady flux level; rods @ 15 
cm; glory hole empty; temp. 
2l°C. 



Steady flux level; 
rods @ 10.1 cm; MRCTC 
set for proper max. 
reactivity conditions . 



0 



Inserted rods @ fast speed. 



Inserted rods @ fast 
speed . 
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PIG. 10. VIEW OF COMPLETED SIMULATOR. WITH 
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FIG. 11. BEAR VIEW OF SIMULATOR CQNSOIE 





0.4 


Rods up; flux permitted to 
increase on 17.5 sec . period 


Rods up; flux per- 
mitted to increase 
on 17.5 sec. period. 


2.1 


Flux level^ 7.6 x 10^ 

(0.76 x 10^ ppa scale read- 
ing); rods retracted to 
20.45 cm; flux permitted to 
decay to 1.0 x 10 . 


Flux level @ 7.6 x 

10^ (76 volts scale 
reading) rods re- 
tracted to 22.5 cm; 
flux permitted to^ 
decay to 1.0 x 10 . 


4„25 


Adjusted rods for steady 
flux level. 


Adjusted rods for 
steady flux level. 


7.50 


Quickly inserted moderator, 
graphite, and lead discs in- 
to glory hole to provide a 
step change of about (+) 

1300 Q of total required 

to make core essentially 
homogeneous , ) 


Quickly rotated MRCTC 
potentiometer to pro- 
vide corresponding 
(+) 1300 py© step change 
in reactivity. 


8.65 


Retracted rods to 19.93 
cm; allowed flux to de- 
cay to 4.0 x 10®. 


Retracted rods to 21.9 cm; 
allowed flux to decay to 
about 4.0 x 10®. 


10.15 


Inserted remainder of 
poly, graphite, and 
lead in glory hole and 
adjusted rods to steady 
flux level of 4.5 x Kr 
f rated power 100 mv. 


Rotated MRCTC potentiometer 
an additional amount corr- 
esponding to about 1300 jy? 
to simulate homogeneous core 
conditions and brought to 
steady flux level of 4.5 x 

10. 


12.05 


Quickly withdrew all 
poly from glory hole to 
provide a step change of 
about ..(-) 2600 ; al-- 

lowed flux level to decay 
to about O.63 x 10®. 


Quickly reset MRCTC to or- 
iginal setting correspond- 
ing to glory hole empty 
conditions thus providing 
a step change in reactiv- 
ity of about (-) 2600 ; 

allowed flux to^decay to 
about 0.54 x 10 . 


13.5 


Inserted rods all the 
way; permitted flux to 
increase to about 4.5 
x 10® and leveled at 
about 4.0 x 10^. 


Inserted rods all the way; 
permitted flux to increase 
to about 4.5 x 10® and ^ 
leveled at about 4.0 x 10 . 



As indicated in Table 3, the position to which simulator rods are retracted 
in each case is different from that given by the rod position indicator on 
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FIG. IE. FLUX-TIME C0KIPARI50IM CURVES 



the reactor console. This implies that the reactivity potentiometer con- 
tributes more reactivity per centimeter of movement than do the actual rods , 
By a trial and error procedure of slightly altering the value of the 
paralleling resistance of the reactivity potentiometer, the actual and 
simulator rod positions might be more closely brought into agreement. 

The sequence of events describing Figure 13 is given in Table 4. 



TABLE 4. SEQUENCE OF EVENTS FOR FIGURE 13 

TIME IN REACTOR SEQUENCE CORRESPONDING SIM- 

M I NOTES ULATOR SEQUENCE 



0 



0.4 



2.1 



3o0 



Approximately steady 
flux level; rods @ 15 cm; 
glory hole empty; temp. 
21°C. 



Inserted rods at fast 
speed. 

Rods up; flux permitted 
to increase on 17.5 
sec. period. 

Retracted rods and 
brought to steady flux 

level of 7.0 x 10 6 . 

Scrammed reactor and 
permitted flux to decay 
to observe effects of 
delayed neutrons. 



Steady flux level; 
rods @ 10.1 cm; MRCTC 
potentiometer set for 
proper maximum reactivity 
conditions . 

Inserted rods at fast 
speed. 

Rods up; flux permitted 
to increase on 17.5 
sec. period. 

Retracted rods and brought 
to steady flux level of 

7.0 x 10 6 . 

Scrammed simulator and 
permitted flux to decay 
to observe effects of 
delayed neutron networks. 



Simulation on scramming the reactor is difficult to reproduce because of 
dashpot action on the rods as they are ejected from the core. This action 
was approximated by manually rotating the reactivity potentiometer to 
the all-rods-out position. As a result, the simulator decay curve on 
scramming appears to be srteeper than the reactor curve. However, the 
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FIG. 13. FLUX-TIME COMPARISON CURVES 



general nature of the two decay schemes is quite similar with the flux 
levels of reactor and simulator being 3.4 x 10^ 2.3 x 10^ respectively 

four minutes after scramming. 

The sequence of events describing Figure 14 is given in Table 5. 



TABLE 5. SEQUENCE OF EVENTS FOR FIGURE 14 



TINE IN REACTOR SEQUENCE CORRESPONDING SIMULATOR 

MINUTES SEQUENCE 



Approximately steady 

flux level; rods® 20.0 
cm; glory hole empty; 
temp. 21 °C. 

0 Quickly inserted all 

poly, graphite, and lead 
discs into glory hole to 
provide a step change from 
about (-) 1300 \ip to 
(+) 1300 pyo conditions. 
Permitted flux to increase 
to about 0.73 x 10 . 



Steady flux level rods 
@ 22.0 cm; MRCTC potentio- 
meter set for proper max. 
react, conditions. 

Quickly rotated MRCTC potentio- 
meter to simulate step change 
of about 2600 y-p bringing 
simulator from a subcritical 
condition of about (-) 1300 ^ 
to a supercritical condition 
of about ( + ) 1300 V-f* . 
Permitted flux level to in- 
crease to about 0.73 x 106 . 



2.5 



Quickly removed all poly 
to provide step change from 
a supercritical condition 
of about (+) 1300 to 

a subcritical condition of 
about (-) 1300 0 



Quickly rotated MRCTC potentio- 
meter back to original setting 
corresponding to a step change 
in reactivity from about (+) 
1300 u.p to about (-) 

1300 . 



The current output in pji amps, from an ion chamber feeds the linear 



recorder on the reactor; the appropriate scales are labeled in the three 
figures. Corresponding to these are the voltage scales on the simulator 
which are also labeled in the figures. For example, the 1 x 10“ UPa 
scale on the reactor corresponds to the 100 volt scale on the simulator, 
the 1 x 10 corresponds to the 10 volt scale, the 1 x 10“"^ upa 
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H|ia to the 0.1 volt scale. The 



to the 1 volt scale, and the 1 x lcf 11 

reading in as well as the reading in volts, of course, are proportional 

to a certain flux level. The flux level at the center of the core at a 

power level of 100 nrw is given as 4*5 x 10^ # This corresponds to a 
0.45 x 10~° lH*>a output from the ion chamber or a 45 volt output from 
the simulator; all other values are scaled accordingly. 

7. Conclusions 

As evidenced by the comparison curves of Figures 12 through 14, reason- 
ably accurate simulation of the time dependent behavior of a specific low 
power thermal reactor may be readily accomplished despite the many simplifying 
assumptions made in the development and electrical duplication of the 
governing transient equations. The final simulator with its relatively 
simple control system provides a useful tool with which to observe typical 
reactor responses resulting from the many sequences of operations to which 
a reactor of this type might be subject. 
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APPENDIX 



CIRCUIT ANALYSIS DETERMINATION OF CURRENT Ig 




Assume current directions as indicated. Therefore, 

Wa = *8 

Now, let A = proportional distance of P from mid or zero position. 
Then 



-E - 0 

h> " R + P/2 (1 + A) 



i 

a 



0 - E 

R + P/2 (1 - A) 



Therefore , 



h ~ R + p/2 (1 + a) + R + P/2 (1 - a) 
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or *8 ' E 



= E 



1 1 

R + P/2 (1 - A) “ R ♦ P/2 (1 + A) 



2 A P/2 

£R + P/2 (1 - A)} £r + P/2 (1 ♦ A )} 

E A P 



R2 



+ R| J~~ (l + a) + (1 - aTJ ♦ <1 - t*) 



h = 



A E 



jr * R + | (1 - A 2 ) 



or 



% = 



A E 



R Q (1 ♦ 1> * fe ft - 
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APPENDIX II 



SCALING FACTOR CALCULATIONS 
Referring to scaling factor equation, namely 



■J 



m - 1 "(vta,! 



[■, 



( 22 ) 



let R^ = 10 hfegohms . This is an arbitrary value and is purposely chosen 

large to permit reasonable values of voltage V which is dependent on the 

very short neutron lifetime of this reactor. 

„ -4 

Next, assume a value of neutron lifetime 1 of 10 sec. Since the reactor 
will not be far from critical at any time, 1* and 1§ will be essentially 

equal to each other. Hence, as a first approximation 1 0 will also be 

-4 

assigned a value of 10 sec. 

Now from control rod calibration data, & can be determined. In view 

/ max 

of the fact that normal operating temperatures are about 23° C, 
will be determined at this temperature assuming a homogeneous core (i.e., 
moderator, reflector, and shielding material in the glory hole). The 
operating range of reactivity here is 

*(-) .034775 to (+) .004550 



Therefore , can be taken as .034775* and an imposed unbalance in the 
resistance values of R will prevent any positive reactivity value greater 
than .004550, In other words, when A is (-) ~ I max) 

A relationship must now be established between voltage E and neutron 
population density n. From known data, 
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<p at 100 mw = 4.5 x 10^ neuts/cm 2 /sec 
and v = 1.2? x 10^ *Y T°K cm/sec 
or v = 2.22 x 10^ cm/sec at 23°C. 

Therefore, n = | = ^ ‘ L^IO ^ = 20,25 neuts / cm3 

at rated power level, 100 mw. 

Because the maximum output of the amplifier is 100 volts and because it 
might be desirable to directly correlate flux with observed voltage, 45 
volts is made to correspond to 4.5 x 10^ neutrons per cm 2 per sec. The 
constant of proportionality 3 is then given by 



J = = _1L. 

n 20.25 



2.22 



Next, a suitable value of potentiometer resistance P is selected so that 
after modifications by an external resistance to a center tap connection, 
a reasonable value of total resistance is still available. A 20K pot 
meets this requirement; and after modification as indicated in Figure 5* 
an overall resistance of 13.33K results (13. OK measured). It is to be 
noted that this value is not that which is to be substituted for P in 
equation (22). A roust be determined which will be compatible with 

the derivation of the reactivity potentiometer relationships given in 
section 4. As far as the circuit is concerned, for zero or critical set- 
ting of the potentiometer the resistances on each side must be balanced. 
Therefore, since a (-) .034775 reactivity is required on one side, or in 
terms of resistance 

( ,, 0 34775 )• loir = ( PpsV(13^ = 11 <5 KSt 

[ .034775 + .00455 J U ' 5; {J/ 1 
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This implies a value of 



P eff = 2 x 11,5 KA = 23.0 KSL 



which should then be utilized in (22) in place of P. Now, because the 
actual resistance of P is only 13.0 K, then 10K additional resistance is 
required in loop 2 as indicated in Figure 6, or R 2 = R, + 10K-fiL . Values 
of are readily available and are given in the summary of scaling factor 
values at the end of section 5. Before calculations begin, a value of 
R must also be selected. First, however, since it is on this resistance 
that the individual values in the six Rc networks will depend, investiga- 
tion to ensure reasonable magnitudes for the Rc delay networks is 
necessary. Based upon this investigation, a final value of R = 119K yields 
suitable magnitudes for all six delay components. Therefore, 

R, = 119K and R 2 = 119K + 10 K = 129 Kji . It is also noted here that R& 
5P which yields better than 1^ linearity between^oand A. Calculations 
can now proceed as follows* 



R(1 + R/P ^ = 119 x 10 3 (1 + 0.034775 



2.55 x io 4 .a 



Now, from (22) 





J 
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_ 1 . _ 10 ^ 

0r ° 7 ' ' 2.55 x 10' 

c ? = .00392 [if 



Likewise , 





J 


J 








Pi R i 




or 


R & ♦ 


. 2.55 x 10 4 




Pi 


Pi 




i 


^ X^sec” 1 


R^megs 




lo 


.00025 14 


102 


0 .0007 


2 


.00085 1.61 


30 


0.0207 


3 


.00241 0.455 


10.6 


0.208 


4 


.00213 0.154 


12.0 


0.544 


5 


.00166 0.0315 


15.4 


2.06 


6 


.00025 0.0125 


102 


0.785 


P = 


2|l p t = .00755 







Next, the method for deciding upon a suitable value of source voltage, 
V, as discussed in section 5> is as follows: 



Power ratio x n 



@100mw 



Now, n 



@30 pw 



30 x 10 



-6 



100 x 10 



-3 



x 20.25 



*"3 3 

or n __ = 6. OP x 10 neuts / cur / sec 

@30 yw 



P 



. K .<r - 1 

K eff 



@ safeties up position = (-) .009625 



from control rod calibration data. From this then, 



eff 



.99047 



or (1-K eff ) 



.00953 



Therefore , 



Sol = .579 x 10 



-4 



and S 0 1 0 *= .5P5 x 10 



-4 



J o- L o 



Therefore , 



V = 



3 S (s °(ma x^ lo) 

R (1 + R/P 



max 



or 



V = 



(2.22) *(10 x 10 6 )-(.5P5 x 10~ 4 ) 
(2.55 x 10 4 ) 



V = 51.0 mv 



It is to be noted that, in the case of the simulation of the AGN-201, 
S 0 and S 0 ( mx ) are the same since it is not desired to vary the source 
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strength (voltage); and therefore only the one magnitude determined above 
is needed to represent the effectively constant Ra-Be neutron source in 
the reactor. 

Finally, if 1* is taken as 10 sec, then the value of S 0 is readily 
found to be 

-4 

S 0 = -- 1 — = ,579 neuts/cm^/sec 

10~ 4 
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